Topological Fermi liquids from Coulomb interactions in the doped Honeycomb lattice 
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We get an anomalous Hall metallic state in the Honeycomb lattice with nearest neighbors only 
arising as a spontaneously broken symmetry state from a local nearest neighbor Coulomb interaction 
V . The key ingredient is to enlarge the unit cell to host six atoms that permits Kekule distortions 
and supports self-consistent currents creating non trivial magnetic configurations with total zero 
flux. We find within a variational mean field approach a metallic phase with broken time reversal 
symmetry (T) very close in parameter space to a Pomeranchuk instability. Within the T broken 
region the predominant configuration is an anomalous Hall phase with non zero Hall conductivity, 
a realization of a topological Fermi liquid. A T broken phase with zero Hall conductivity is stable 
in a small region of the parameter space for lower values of V . 
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Introduction. — The quantum Hah effect [l| and the 
anomalous Hall (AH) effect [2] have given birth to a 
new paradigm in condensed matter based on momen- 
tum space topology [3]. Very interesting developments 
have followed based on the recent ideas of getting a 
Hall conductivity or Landau levels without external mag- 
netic fields 0, m. These ideas have given rise to new 
areas of research and associated new materials as the 
topological insulators (6| and the even more interesting 
topological metals 0- These systems allow the realiza- 
tion of beautiful fundamental ideas shared by different 
branches of physics like charge fractionalization or Ma- 
jorana fermions [8, 9]. 

The interplay of the underlying lattice and the elec- 
tronic interactions plays a very important role in the 
physics of these systems. In topological insulators the 
spin-orbit coupling is the main ingredient to get non triv- 
ial topological phases. In the topological metals time re- 
versal symmetry (T) breaking without a magnetic field is 
the key ingredient, which can be realized through current 
(or bond) ordering: the electrons spontaneously form 
current loops, which interact among themselves in such a 
way that the state is self-consistently maintained. These 
phases were discussed in other contexts in lol- 13|. 

Because of its special topology and long before the syn- 
thesis of graphene, the Honeycomb lattice has played a 
predominant role in the modeling of topological states of 
matter. One of the earliest examples of a quantum Hall 
effect without a uniform externally applied magnetic field 
is due to Haldane who obtained a T broken state in a 
tight binding model in the Honeycomb lattice with com- 
plex values of the next to nearest neighbors hopping pa- 
rameters. Ever since, the search for realization of sponta- 



neous (AH) effect has been very intense in the literature. 
The problem turns out to be very hard and the proposed 
models usually involve the inclusion of hopping or inter- 
actions beyond the nearest neighbors as in the original 
Haldane model, or very e labora ted lattice structures as 
the Kagome or pyroclore 
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In this work we show that an AH phase exists as a sta- 
ble ground state in a simple nearest neighbor tight bind- 
ing model with nearest neighbor Coulomb interaction in 
the Honeycomb lattice. The T broken state support- 
ing this phase is made possible by enlarging the original 
two atom unit cell to a six atom unit cell (see Fig. [T]) 
which allows for local non-zero current states and thus T 
broken phases. Our new phases are a realization of the 
topological Fermi liquids described in 0, [HI . 

The model. — The extended Hubbard Hamiltonian that 
we consider for spinless fermions in the honeycomb lattice 
reads 

H = -t^^alhrJrd + V^^alarhlj^^hrJrd + h.c. , (1) 



Y,8 



r,(5 



where t is the nearest neighbor hopping and V the near- 
est neighbor Coulomb repulsion. We use standard no- 
tation where ar (hr) annihilates an electron at position 
r in sublattice A (B). The two inequivalent sublattices 
A and B are depicted in Fig. [U^a), along with the ba- 
sis vectors ai = ^(-1,^3) and as = ^(l,\/3) for 
the case of a two atom unit cell. The vectors 5 refer 
to the three vectors connecting nearest neighbor sites, as 
shown in Fig. [IJa). To allow for T broken phases as mean 
field solutions we use an enlarged unit cell, containing 
six atoms, which also permits Kekule type distortion as 
shown in Fig.fljb). The basis vectors of the enlarged cell 




FIG. 1: (color online), (a) Two- atom unit cell and example 
of uniaxial distortion, (b) Six atom unit cell and Kekule dis- 
tortion, allowed in the enlarged unit cell, (c) A pictorial rep- 
resentation of the nine complex order parameters considered 
in this work in the mean field decoupling of the Hamiltonian 



in real space are ai = ^(— \/3, 1) and a.2 = ^(\/3, 1), 
and the respective unit cell vectors in reciprocal space 
= 5afe(-l'V^) h, = 5^(1, V3). With 
this choice the unit cell in direct space is three times 
bigger while in reciprocal space the Brillouin zone (BZ) 
becomes folded, i.e. it is three times smaller. This gives 
rise to a tight binding model whose wave function is a six 
component spinor of the form 
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[al (k), b\(k), 4(k), bl(k), 4(k), 4(k)]. (2) 



Since we are interested in the electronic phases with 
broken T we do not consider for the time being charge 
ordered phases. Under these conditions the most gen- 
eral mean field Hamiltonian depends on nine complex 
parameters which can be grouped in a 3 x 3 matrix, 
and that can be shown to be /c— independent. The mean 
field equations can be written in terms of the mean field 
averages of the form (bU'k) a i {]<.)) mf as 



7L^'{6](k)a,(k))MF, 



(3) 



where is the number of unit cells, 7k is a 3 x 3 matrix 
given by 
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and the momentum sum runs over the folded BZ. The 
nine complex order parameters of our mean field de- 
coupling represent the nine bonds in the enlarged unit 
cell, as pictorially represented in Fig. Hl^c). We solve 
Eq. (|3]) self-consistently with the constrain imposed by 
the Luttinger theorem [2o|, which reads (ignoring loga- 
rithmic corrections in fermion number Ne)^ 



(4) 



kJ 



where n is the electron density per unit cell relative to 
half filling (which in our case corresponds to n = 0), 




FIG. 2: (color online). Mean field phase diagram. Legend: 
(S) symmetric phase, i.e. bare graphene with a uniform renor- 
malization of the hopping; (K) Kekule distortion with hop- 
ping renormalization as shown in the inset; (P) Pomeranchuk 
distortion of the Fermi surface and hopping renormalization 
as shown in the inset; (K+P) coexistence of Kekule and 
Pomeranchuk distortions; (T-I) and (T-II) anomalous Hall 
phases discussed at length in the text; (RS) broken symmetry 
state with real hopping parameters, the distortion is neither 
Kekule type nor Pomeranchuk (reduced symmetry). 



£i(k) is the mean field dispersion for the / band, and 
nF[£|(k), /i] is the Fermi distribution function. From 
Eq. (jl]) we get the renormalized chemical potential /i self- 
consistently. 

The phase diagram and the AH phase. — The mean 
field phase diagram is shown in Fig. [2] where we plot the 
different phases (defined in the caption) as a function of 
the interaction strength V in units of the hopping pa- 
rameter t and the electron density n. The density varies 
from n = (half filling) to well above the VH filling 
which with our convention occurs at nyu = 0.75 . At 
each point in the phase diagram a mean field Hamiltonian 
can be extracted which can be seen as a free Hamiltonian 
with new effective hopping parameters renormalized by 
the interaction. At low values of V the symmetric phase 
(S) represents standard graphene with a uniform renor- 
malization of the hopping. Close to half filling for in- 
creasing values of V slightly above y = 2t we recover 
the Kekule phase (K) described in [sl which evolves to 
a Pomeranchuk phase (P) through a finite coexistence 
region (K+P). Our calculation shows that a standard 
Pomeranchuk instability with an anisotropic renormal- 
ization of the hoppings as shown in the inset, is a very 
robust phase around the VH filling from zero to high 
values of V. This result was already obtained without 
the BZ folding in [2l|, |22|. The preferred phase is a ne- 
matic one where the Cq symmetry of the original lattice 
is broken to a C2 . The inset shows one of the three equiv- 
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FIG. 3: (color online). Pictorial representation of the or- 
der parameters corresponding to the T broken phases T-I (a) 
and T-II (b) discussed in the text. The thickness of the bonds 
represent the modulus of the hopping parameter and the di- 
rection of the arrows represents the sign of the phase when it 
has a complex value. A bond without an arrow means a real 
hopping. 




FIG. 4: (color online). (Left) Two lowest energy bands for 
the mean field Hamiltonian in the AH phase obtained with 
V — 5t, n = 1.13. The hexagonal line marks the position of 
the Fermi level. (Right) The Fermi surface. 



alent configurations oriented along the crystal principal 
directions. The phase named reduced symmetry (RS) oc- 
curring at higher values of the electron density and the 
interaction is a broken symmetry state with real hopping 
parameters. The distortion is neither Kekule type nor 
Pomeranchuk. 

The novel topological Fermi liquid phases appear near 
n = 1. There are two T broken phases labelled T-I and 
T-II in Fig. [2] which are the most stable configurations 
just above the VH filling for moderate values of V be- 
ginning dit V ^ 3t. They are pictorially described in 
Fig. ISl^a) and[3l^b), respectively, where the nine complex 
order parameters of our mean field decoupling are shown. 
The direction of the arrows represents the sign of the 
phase of the given complex hopping, and the thickness of 
the line represents its modulus. The phases can also be 
understood as patterns of orbital currents. Current con- 
servation at each of the six atoms in the unit cell plus the 
zero overall fiux condition allow for only two independent 
T breaking phases, T-I and T-II in our notation, defined 
by their corresponding flux pattern in the unit cell. As it 
is clear from Fig. 3 both phases are only possible if the 
unit cell is enlarged. We note that in addition of having 
the non-trivial fluxes described their structure includes a 
Kekule distortion of the bonds. 

The discrete symmetries of the mean field Hamiltonian 
help to classify the topological properties of a given phase 
|l9| . Following [23| the discrete symmetry operations T 
and inversion symmetry (X) amount to the transforma- 
tion, 

with 64 = hi. The phase T-I breaks T and X, but pre- 
serves TX. The T-II breaks T but X is preserved. 

In what follows we will show that the T-II phase that 
dominates the T broken part of the phase diagram is 
indeed an AH phase. As discussed in the (non- 
quantized) Hall conductivity of the topological metals 
as the ones encountered in this work is a property of the 



Fermi surface. At a given point of the phase diagram 
it can be computed from the single particle Bloch states 
|^/(k)) associated to the appropriate mean field Hamil- 
tonian from the expression: 

= ^^^l^f (k)n^h(k),M], (6) 

k,l 

where V is the volume of the unit cell, and l]^^(k) is 
the Berry curvature defined from the Berry connection: 

^r(k) = -z(^Kk)|V^i,)^Kk)), ^f(y) = VMf(k) - 

V^^f (k). The T-II phase is of the type II in the classifi- 
cation given in [19]: it breaks T but preserves I and the 
Hall conductivity is generically non zero. The T-I phase. 
Fig. Ilfa), breaks T and X but preserves TX so it corre- 
sponds to a T broken phase of type I and has zero Hall 
conductivity. We have further confirmed this picture by 
numerical computation of the Hall conductivity Eq. ([6]). 

A very neat analysis of the topological properties of 
the various metallic phases in the phase diagram can be 
done by studying the low energy effective bands. We 
have plotted the mean field band structure of the T- 
II phase in Fig. |4] obtained for the parameter values 
V = 5t, n = 1.13. Focusing on the relevant bands around 
the Fermi level, it is easy to understand the qualitative 
behavior of the non vanishing Hall conductivity in this 
phase. The Fermi level crosses a massive Dirac structure 
around the F point and so there is a non-zero contri- 
bution to the AH conductivity of this cone. The non- 
quantized contribution to the AH conductivity from the 
cone is given by i25i] 

^ Min,V) 
'"^ 2h\tl{n,V)\ 

where M(n, V) is the gap at the F point and /i(n, V) is 
the renormalized chemical potential relative to the mid- 
dle of the gap, both of which depend strongly on the 
parameters of the phase diagram. 
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To better understand the nature of the T broken 
phases we note that they arise in the region of the pa- 
rameter space close to the density where there are four 
electrons per unit cell: n = 1. This is a very special fill- 
ing: Not only it is commensurate with the lattice, but it 
enhances the formation of current loops self-consistently 
maintained at each hexagon following the configurations 
shown in Fig. [3l We have seen that in the T broken part 
of the phase diagram, along the line n = 1 the system be- 
comes an insulator. The band structure near the Fermi 
level is similar to the one shown in Fig. |4] left but the 
cones are further apart and the Fermi level lies in the 
gap. Away from this line we have the situation described 
before. The majority of electrons will still form currents 
as these in Fig. [3] and the excess (defect) electrons are re- 
sponsible for the metallicity of the system. This picture 
allows also to understand the asymmetry of the phase 
diagram with respect to the line n = 1 (see Fig. [2]) as a 
result of the competition between the soft Fermi surface 
associated to the high density of states near the VH fill- 
ing and the orbital ordering driven by high to moderate 
values of the interaction V at the commensurate filling. 
Since the density of states (DOS) is very steep around 
this density, moving to the left (for lower values of n) 
the DOS grows very fast towards the VH filling and the 
Pomeranchuk instability takes over. Going to densities 
n > 1 the DOS decreases and the local currents respon- 
sible for the T broken phases survive in a larger region. 

Discussion and future. — Part of the physics discussed 
in this work can be tested in actual graphene samples. 
The simple deformation of the Fermi surface pointing to 
a Pomeranchuk instability is a very robust phase that 
may prevail even if other instabilities not considered in 
this work are allowed. The AH phase can be more diffi- 
cult to observe in graphene since it occurs at higher val- 
ues of the interactions and is accompanied by a Kekule 
distortion but it could potentially be tested in cold atom 
experiments with optical lattices [26|, [l^ . 

Other phases may compete with the ones described in 
this work when charge decoupling and spin degrees of 
freedom are included in the system. It will be interesting 
to see how they compete with the AH phase obtained in 
this work. Very appealing possibilities will open in the 
Pomeranchuk region of the phase diagram in Fig. [2] when 
spin is included along the lines of [28|. Spin effects have 
also been explored recently in [29|. A detailed study of 
the physical properties of the various phases encountered 
as well as the detailed nature of the phase transitions will 
be examined in a forthcoming publication. 

Conclusions. — We have found a spontaneous symme- 
try breaking to an AH phase in a tight binding model in 
the Honeycomb lattice with only nearest neighbor hop- 
ping parameters and Coulomb interaction. The extra 
physics required to get such a phase is provided by the 
folding of the BZ that allows for spontaneous non-zero 
currents with zero overall magnetic flux to form inside 



the unit cell generating T broken phases. The T broken 
phase is predominantly an AH metal of the type H in 
the classification given in [l9| where the interaction V 
gives rise to orbital current fluctuations together with a 
Kekule distortion. 

The findings of this work open a whole set of possibil- 
ities for new realization of exotic phases based on lattice 
models. Enlarging the unit cell is a very simple proce- 
dure that increases enormously the phase space of any 
given lattice. This is exemplified in the model studied 
here where in addition to the AH phase we have found a 
very rich phase diagram even when neglecting spin and 
charge instabilities. 

This research was supported by the MEG (Spain) 
through grants FIS2008-00124, PIB2010BZ-00512. F.J. 
acknowledges support from NSF through Grant No. 
DMR-1005035. 



[1] Z. F. Ezawa, Quantum Hall Effects - Field Theoretical 
Approach and Related Topics (World Scientific, Singa- 
pore, 2008). 

[2] N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and 
N. P. Ong, Rev. Mod. Phys. 82, 1539 (2010). 

[3] G. E. Volovik, The universe in a helium droplet (Claren- 
don Press, Oxford, 2003). 

[4] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

[5] F. Guinea, M. I. Katsnelson, and A. G. Geim, Nature 
Physics 6, 30 (2010). 

[6] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

[7] F. D. M. Haldane, Phys. Rev. Lett. 93, 206602 (2004). 
[8] C.-Y. Hou, C. Chamon, and C. Mudry, Phys. Rev. Lett. 

98, 186809 (2007). 
[9] M. Franz, Physics 3, 24 (2010). 
[10] A. A. Nersesyan, Phys. Lett. A 153, 49 (1991). 
[11] E. H. Lieb, Phys. Rev. Lett. 73, 2158 (1994). 
[12] S. Chakravarty, R. B. Laughlin, D. K. Morr, and 

C. Nayak, Phys. Rev. B 63, 094503 (2001). 
[13] C. M. Varma, Phys. Rev. B 73, 155113 (2006). 
[14] K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Phys. 

Rev. Lett. 103, 046811 (2009). 
[15] S. Raghu, X.-L. Qi, C. Honerkamp, and S.-C. Zhang, 

Phys. Rev. Lett. 100, 156401 (2008). 
[16] T. Pereg-Barnea and G. Refael, arXiv: 101 1.5243 (2010). 
[17] J. Wen, A. Ruegg, C.-C. J. Wang, and G. A. Fiete, Phys. 

Rev. B 82, 075125 (2010). 
[18] Z. Qiao et al.. Physical Review B 82, 161414(R) (2010). 
[19] K. Sun and E. Fradkin, Phys. Rev. B 78, 245122 (2008). 
[20] J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 
[21] B. Valenzuela and M. A. H. Vozmediano, New J. Phys. 

10, 113009 (2008). 
[22] C. Lamas, D. Cabra, and N. Grandi, Phys. Rev. B 80, 

075108 (2009). 

[23] J. Manes, F. Guinea, and M. A. H. Vozmediano, Phys. 

Rev. B 75, 155424 (2007). 
[24] C. Kane and E. Mele, Phys, Rev. Lett. 95, 146802 (2005). 
[25] A. Cortijo, A. G. Grushin, and M. A. H. Vozmediano, 

Phys. Rev. B 82, 195438 (2010). 



5 



[26] L. B. Shao, S.-L. Zhu, L. Sheng, D. Y. Xing, and Z. D. [28] C. Wu, K. Sun, E. Fradkin, and S. C. Zhang, Phys. Rev. 

Wang, Phys. Rev. Lett. 101, 246810 (2008). B 75, 115103 (2007). 

[27] K. L. Lee, B. Gremaud, R-Han, B.-G. Englert, and [29] T. Li, arXiv: 1103.2420 (2011). 

C. Miniatura, Phys. Rev. A 80, 043411 (2009). 



